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We give two formulations of exclusion statistics (ES) using 
a variable number of bosonic or fermionic single-particle states 
which depend on the number of particles in the system. As- 
sociated bosonic and fermionic ES parameters are introduced 
and are discussed for FQHE quasiparticles, anyons in the low- 
est Landau level and for the Calogero-Sutherland model. In 
the latter case, only one family of solutions is emphasized 
to be sufficient to recover ES; appropriate families are speci- 
fied for a number of formulations of the Calogero-Sutherland 
model. We extend the picture of variable number of single- 
particle states to generalized ideal gases with statistical in- 
teraction between particles of different momenta. Integral 
equations are derived which determine the momentum distri- 
bution for single-particle states and distribution of particles 
over the single-particle states in the thermal equilibrium. 

PACS Numbers: 05.30.-d, 05.70.Ce, 74.20. Kk, 73.40.Hm 



The idea of reduction of single-particle Hilbert space 
upon adding particles into the system introduced by HaL 
dane into definition of statistics for identical particlesu 
has recently received great attention. This idea was orig- 
inally proposed to give an alternative definition of frac- 
tional statistics based on a generalized exclusion principle 
[exclusion statistics (ES)], with a specific law of reduc- 
tion of the Hilbert space determined by an ES parame- 
ter. Associated formula for the dimension of the many- 
particle Hilbert space was suggested!!! which was then , 
plied to FQHE auasiparticles to evaluate numerically, 
and analyticallyta their ES parameters in simple cases. 

The latter formula applied locally in the phase space 
was also used as a starting point for formulation of the 
statistical mechanicscm (see also Ref. 7). The statisti- 
cal distribution fop, ES turned out to coincide with that 
previously derived^ for fractional statistics in one dimen- 
sion [the latter was introduced in the algebraic (Heisen- 
berg) approach to quantization of identical particlesu and 
was arguecffij to be modeled by systems with an inverse 
square interaction (the Calogero-Sutherland modelsO)]. 
On related studies of ES in terms of the Calogero- 
Sutherland models, see Refs. 12-15. 

Interpretation of anyons confined to the lowest Lan- 
dau level (LLLV,i«, terms of ES statistics was discovered 
in various waya9ll3~£3. This is consistent with the inter- 
pretation of anyons in the LLL, where their dynamics be- 
comes effectively one-dimensional, in terms of fractional 
statistics in one dimension in the algebraic approach!!!!. 

The idea of reduction of the space of single-particle 



states was used to define statistical interaction between 
particles of different momentatj. The latter notion was 
also discussed in the language of generalized single=state 
statistical distributions for free identical particlesEl 

Haldane has defined ES by the condition that upon 
adding AN particles into the system (we restrict our- 
selves to the case of one species of particles), the number 
of single-particle states D available for further particles 
is altered by 



AD = —g AN , 



(1) 



where g is the ES parameter (with g = and g = 1 for 
bosons and fcrmions). He also suggested the formula for 
the dimension of the space of many-particle states: 



W = 



(D + N - 1)! 
N\(D-1)\ 



(2) 



If one denotes by G the number of single-particle states 
available for the first particle and uses the expression 



D = G — g(N - 1) 



(3) 



consistent with (|l|), then (|^) interpolates between the 
numbers of ways of placing TV bosons (g — 0) and N 
fermions (g — 1) over G single-particle states. 

We stress that in order to derive the dimension of 
the space of many-particle states starting from the di- 
mension of the space of single-particle states, it is nec- 
essary to know "accommodation" properties of single- 
particle states. In the usual interpretation of Eqs. (0)- 
@ as counting ways of placing N particles into G states, 
Eq. ^ (which in fact determines implicitly "accommo- 
dation" properties of single-particle states) should be 
viewed as independent of ([!]). 

In this paper we give an alternative interpretation of 
Eqs. (|l|)-(||), in terms of usual bosonic single-particle 
states. In this case, @ is a consequence of (Q). We 
also give a dual description of ES which uses fermionic 
single-particle states. The two formulations of the ES are 
somewhat similar to the boson-based and fermion-based 
descriptions of anyons. We extend these formulations 
to involve local statistical interaction in the phase space 
as well as statistical interaction between particles of dif- 
ferent energies. Based on this, we reexamine systems 
of anyons in the LLL, the Calogero-Sutherland systems 
and interacting systems of spinless particle solvable by 
the thermodynamic BetheL-apsa^JlTBA), refining and 
extending previous resultsaEHljTJ'tj on thermodynamic 
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equivalence of the above systems to systems of free iden- 
tical particles. 

Bosonic and fermionic single-particle states. — Con- 
sider a system that may be described in terms of bosonic 
(fermionic) single-particle states and assume that the to- 
tal number of single-particle states depends on the num- 
ber of particle in the system: D b ^ = D b ^{N). Single- 
particle states may thus appear or disappear on adding 
particles into the system. The formula similar to ([j]) 
then defines variation the number of single-particle states 
upon adding AN particles into the system: 



AD b - f (N) = -g bJ AN 



(4) 



where we have introduced the bosonic and fermionic ES 
parameters g b and gf . In place of (0), we have 



D b ' f = G-g b ' f (N- 1), 



(5) 



where G is the number of single-particle states in the 
absence of particles. The dimension of the space of many- 
particle states follows from to be 



W 



{D b + N - 1)! 
N\(D b -l)\ 



W f 



N\(Df - N)\ 



(6) 



Note that W b and W s with D b and D f from (|) coincide, 
resulting in the same statistical properties, if 
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fJ J 



(7) 



We see that Haldane's original ES parameter corre- 
sponds to the bosonic ES parameter in our notation. 

As an example, we consider ES parameters for FQHE 
quasiparticles. The ES parameters for the quasiholes 
(qh) and quasielectrons (qe) above the Laughlin states 
with the filling factors v = — with m odd were found by 
estimations of the dimension of many-particle space for 
quasiparticles in the spherical geometryn ES param- 
eters were identified by comparison with (Q)-^), which 
corresponds to the bosonic ES parameters. The results 
of Refs. 2-4 then read g b qh = ± and g b qc = 2 - Ac- 
cording to (j?]) , we may as well say that the fermionic ES 
parameters for the qh and qe are 



5qo 



= "3& = 1 



1/m 



possessing the symmetry g^ e = —g qh which was originally 
suggested for the bosonic ES parameters!! 

We thus have the two descriptions of ES, in terms of 
bosonic and fermionic single-particle states, respectively. 
This is reminiscent of the two descriptions of anyons in 
two spatial dimensions: anyons with the same exchange 
statistics parameter 9 with < 9 < 2tt can be described 
either as bosons carrying the magnetic flux (in units of 
quantum flux 0o = 27r/e, e > 0) 4> b = 9/tt with < 
cf) b < 2 or as fermions carrying the magnetic flux <jy with 
— 1 < cf>f < 1, where the two descriptions are connected 
by ~ 



4> b = l + 4> f . (8) 

This relation is quite similar to (0) . We show below that 
for anyons confined to the LLL where the notion of ES ap- 
plies, this similarity becomes the precise correspondence: 
the bosonic and fermionic ES parameters for anyons in 
the LLL exactly coincide with <j) b and $ . 

Particles of the same energy and anyons in the LLL. 
— Let all the particles have the same energy e. It follows 



from (|6|) that the entropy S — In W iscli 

S bJ = D b ' f [±(l±n b ' f )la(l±n bJ ) -n b ' f hin b ' f ] (9) 

(the upper and lower signs refer to the bosonic-state and 
fermionic-state descriptions, respectively), where 



i bJ ee N/D bJ = n/(l - g bJ n) , 



(10) 



and n — N/G. We maximize with respect to n b 'f 
subject to the constraints of fixed total number of parti- 
cles N = D b >fn bJ and the total energy E = D b < f sn b ' f , 
with the associated Lagrange multipliers /3/i and —f3 
{(3 = 1/T, and /j, is the chemical potential). With (|lC|), 
we obtain the equation for n b ^ in the thermal equilib- 



n b ^(l±n bJ ) ±9 



(11) 



where x — e^^~ E \ 

The equation of state can be derived from (^) and ([l]^ 
using the relations fi = E — TS — fiN and fl = —PA 
(specified for two spatial dimensions with A the area oc- 
cupied by the gas) to be 



P/3 = ±(G/A) ln(l±n fc ' / ) 



(12) 



Let us compare this equation of state with—that ob- 
tained for anyons in a strong magnetic fieldEil, where 
only particles in the LLL (all having the same energy) 
contribute into the equation of state. Regarding anyons 
as bosons carrying the magnetic flux of value 4> b (0 < 
cf) b < 2) in the direction antiparallel to the external mafc 
netic field B, the equation of state was found to beEj 
Pj3 = p L ln[l + v/(l - 4> b v)], where p L = eB/2n is the 
density of states in the LLL, v — N/p^A is the filling 
factor. We may as well rewrite this equation of state 
regarding anyons as fermions carrying magnetic flux <j>* 
(—1 < 4>f < 1) connected with <j) b by (||). The two equa- 
tions of state coincide with (|l^), if one identifies the den- 
sity of states G/A in @ with pl (and hence, n with v), 
and, in addition, makes the identifications 



h b <f 



(13) 



Let us comment on formula ^ in this context. Con- 
sider the mean-field approximation corresponding to the 
change of a flux carried by a particle with a uniform mag- 
netic field of the same flux. Adding a particle into the 
system then diminishes the total magnetic flux by 4> b ^ , 
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hence the the number of states in the LLL, equal to the 
total magnetic flux of in units of the flux quantum, is 



D bJ = eBA/2ir - 4> bJ N . 



(14) 



This corresponds to (|^) with g b 'f from (|l3|). But we stress 
that (14) holds only in the mean-field approximation. 

Statistical interaction between particles of the same 
momentum. - Let now particles in different single- 
particle states may have different energies. For definite- 
ness, we assume that single-particle states may be labeled 
with the momenta which particles have in these states. 
We apply the construction of the previous section to par- 
ticles of the same momentum fc (supplying all the exten- 
sive quantities of the previous section with the subscript 
fc, which turns all the quantities into densities per unit 
momentum). We assume that upon adding particles of 
momentum fc into the system, there may appear or disap- 
pear single-particle states of only the same momentum: 



D b k J = G k - g bJ N k 
The entropy for the whole system is 



S bJ = Y,G k InW; 



k 



(15) 



(16) 



where are obtained from ([|) by the changes N — > 

N k and D b 'f -> D b k ' f . A gain, the statistical properties 
given by the bosonic and fermionic pictures are identical 
if the condition (Q) is fulfilled. 

For riil (= N)./D, ) in equilibrium, in place of (11), 
we have 



Xk , 



(17) 



where Xk — £k \ and Sk is the energy of a particle 

= n k /(l - 



with momentum fc. The substitution 



~~i J — - 



g b Tik), yields the equation for nk (= Nk/Gk)'- 



n k (l - g bJ n k )^ bJ ' (I - g bJ n k ±n k ) 



±9 



i>,/_ 



x k - (18) 



For the bosonic ES parameter, ( |18[ ) recovers the equa- 
tion that was derived in the formulation of ES where 
(||) was viewed as |eounting the ways of placing N par- 
ticles over G stateaS'El. That formulation implied invari- 
able number of single-particle states, and the function nk 
played the part of the distribution of particles over states 
fc ('statistical distribution'). Unlike that, in the present 
formulation of ES, with variable number of single-particle 
states, the distribution of particles over states is given by 
the functions n k . 

Eq. ([Tsl ) for the bosonic ES parameter was discussed in 
the context of thejQalogero-Sutherland model as a system 
which reveals ESHOlij. Here we extend that analysis 
with specific emphasis on the fact that only one family 
of solutions of this model is sufficient to recover ES. 

We discuss the Calogero-Sutherland model in a har- 
monic well governed by the Hamiltonian 



H = - 



2 f-' dx 2 

l — l 1 



A(A-l) 



-UJ' 



(19) 



Consider the coordinate domain x\ < • • • < xn- Due to 
the symmetry of the Hamiltonian ( |l9| ) under the change 
A — * 1 — A, the system (|l9|) has two classes of solutions 
of the form 

V/(A) = A A <MM, A), ^ n (A) = A^Mta}, 1 - A), 

(20) 

where A = ni>j( x i — x j)> i x i\ 1S the set with I = 
1,...,N, and *+({*,}, A) - P({%i}, A) exp(-| £\ x i)> 
where P({xi},\) is a polynomial symmetric in the par- 
ticle coordinates {x{\ with a non-zero free term. The 
solutions ip (X) and ip ll (X) are quadratically integrable 
and non-singular for A > and for A < 1, respectively. 

We may rewrite the solutions of class II as ip u (X) = 
A- A $_({x/}, -A), where -A) = A<P+({xi}, 1- 

A) is an antisymmetric function of the coordinates. We 
can then extend continuously the solutions tp 1 for A > 
and ifj U for A < 1 to the regions xq\ < • •• < xqn, where 
Q is a permutation of 1, . . . , N, in the symmetric and an- 
tisymmetric way, respectively. Making, in addition, the 
changes A — * Xb in t/j 1 and — A — * \p in ip u , we arrive at 
symmetric ('boson') and antisymmetric ('fermion') wave 
functions tp B an d ip F °f the- form 



B,F 



= IAI 



*&±({xi},Xb,f)-. 



(21) 



with As > and A^ > — 1. These function are solu- 
tions of the boson and fermion Hamiltonians H B and 
H F which are obtained from the Hamiltonian ( |l9| ) by 
the changes of the two-body potential V{x) with the po- 
tentials 



V*'*{x) = \b,f(\b,fTI)/x 2 



(22) 



where Xb > and Xp > — 1. 

The energy levels for the solutions ip 1 and V 11 or ; 
equivalently, the solutions ip B and ip F (21), are given 



by E B ' F = E B ' F + J2l=i ^ n U where {n{\ is a set of 
non-negative integers obeying < m < ■ ■ ■ < njv, 
and E B = oj[\N + ±X B N{N - 1)] and E F = u[\N + 
|(1 + Xf)N(N — 1)] are the ground state energies (see 
also Rcf. 23). The corresponding partition functions are 
(q = e-^) 



N 



7 B,F 



IK 1 

i=i 



(23) 



The family of solutions tp B ( f2l| ) for the boson Calogero- 
Sutherland model (or, equivalently, the family of solu- 
tions of class I for ( |l9|) ) was analyzed in Ref . 8 in the limit 
lo — > in terms of generalized single-state partition func- 
tions and statistical distributions in the approach with 
nonvariable number of single-particle states, resulting in 
the statistical distributions (f8h with g b = Xb- 
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Repeating the analysis of Ref. 8 for the family of so- 
lutions tp F ( |2l| ) for the fermion Hamiltonian (or, equiv- 
alently, for the solutions of class II in (20)), with the 



(30) 



partition functions Zjy (j23|), we arrive at the statistical 
distributions n k coinciding with those given by ( |l8| ) with 
9 f = Af- 

We note that the analysis given in Ref. 13 implies using 
the family of the fermion solutions tp F of (|2l]) with the 
parameter g used in Ref. 13 equal to our Xp + 1. Refs. 12 
and 15 use the family ip B (|2l|). 

Statistical interaction between particles of different mo- 
menta. — We now give a generalization of the consid- 
erations of the previous section assuming that on adding 
particles of momentum k into the system, there may ap- 
pear or disappear single-particle states of any other mo- 
menta (cf. Eq. (ph): 



D k f = G k -J2Gk'9 bJ {k,k')N k , 



(24) 



where g b ''(k, k') may be called the bosonic and fermionic 
statistical interaction functions. 
The entropy (|g|) is written as 

S bJ = J2 Gk[±(d b k J ± n k ) ln(d b k J ± n k ) 

k 

Td h k J lnd b k J -n k lnn k ], (25) 

where 

d b k f = D b k f /G k = 1 - £ G k ,g bJ {k, k')n k , . (26) 

k' 

It follows from (p5|)-(p6|) that the bosonic and ferpionic 
pictures result in the same statistical properties i£j 



2tt 



(k,k') = —5(k-k')+gf(k,k') 



(27) 



generalizing (^|). 

Maximization of the entropy with respect to n k , sub- 
ject to the total particle number and total energy con- 
straints, yields the equation for n k in equilibrium 

ln^^J2 G '"{ ± 9 bJ (k',k)lnd b k / 

Xk ft' 

+ [^5{k' -k) T 9 bJ (k', k)] ln(4/ ± n k ,)} , (28) 

with d b k ^ of (^6|). Eq. (^8|) can be rewritten as an equa- 
tion for the distribution of particles over single-particle 



states n 



bj. 



In = £ G k , [^5(k> -k) T 9 bJ (k', k)] HI ± n b k /). 

k' 

(29) 



which follows from ( p6|) and the equality n k = d b k n k , 
enable one in principle to determine all the ther- 
modynamic quantities of the ideal gas (e. g. E = 
E fe G k d b k f e k n b k J ). For g b >f(k, k 1 ) oc 6{k - k'), Eqs. (||) 
and (^) reduce to Eqs. ( |l8| ) and ([l7]), respectively. 

The fermionic state description is convenient to study 
low-temperature properties since n k has a non-singular 
simple form for T = 0. Indeed, for T = 0, all the states 
with |fc| < fco, where ko is some boundary momentum 
(the 'Fermi momentum'), are occupied, the others are 
empty, that is, n k = 1 for \k\ < fc , and n k — for 
|fc| > ko- From Eq. ( |30| ) we then get for T = the integral 
equation for the distribution of the fermionic states 

d{ = \- J2 G k ,gf(k,k')d(,, (31) 
fc'|<ft 

as well as the equation for n k ( — ko < k < ko) 

n k = l- G k ,g f {k,k')n kl . (32) 

fc'|<fto 

f 

For sufficiently low temperatures, n k will slightly differ 
from n{| T =o- 

Refs. 19 and 14 argued that the TBA equations for 
spinlcss particles may be viewed as encoding statistical 
interaction between particles of different momenta. We 
now discuss relation of the above single-particle state def- 
inition of statistics to state counting in the TBA. 

In the description of the TBA we follow Refs. 25 and 
26. Let the wave function have the Bethe ansatz form 



V(xi,. . .,x N ) 



N 

A(P) exp (i kpj Xj ) 



p 



3=1 



(33) 



in the asymptotic region x\ <C ■■■ <C xn, where 
ki > ■■■ > fcjv, and P is a permutation of 1,2... A. 
The coefficients A(P) = A(kpi, . . . , fcpjv) are related 
by A(...k',k...)/A(...k,k' ...) = S(k - k'), where 
S(k — k') = — e - lf) ( k - k ) i s the two-body 'scattering ma- 
trix', and 0(k — k') is the two-body 'phase shift', an- 
tisymmetric in k — k' , with k and k' the particle mo- 
menta. Imposing periodic boundary conditions on the 
wave function (|3^) yields the Bethe ansatz equations 
gjftji fjj.. S(kj — k{) = 1 for allowed momentum val- 
ues. Note that the wave function (B3l), which vanishes 



^ /Tttk . • , i r? — p + ph) over which the particle; 

Eq. (|2S)) together with the equation for the distribution determined by the integral relation i 



if kj = ki for any two momenta (5(0) = — 1), implies 
description of the system in terms of fermionic states. 

In the thermodynamic limit the Bethe ansatz equa- 
tions results in the following picture. The density of lev- 
els rf (being the sum of the particle and hole densities, 
Ph) over which the particles are distributed is 



of single-particle states d' 



k ' 
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2nr f (k) = 1 - J y{k- k')p{k')L dk 1 , 



(34) 



with ip(k) = j^dO{k) / dk the derivative of the phase shift, 
thus depending on the particle density. Fermionic nature 
of the states is also displayed in the expression for the 
entropy: 



-(r^ — p) ln(r^ — p) lnr^ — p\np\Ldk. 



We see that the TBA picture exactly corresponds to 
the above definition of statistics in terms of fermionic 
states: Eq. ( |34| ) and the entropy coincide with (^6|) and 
( p5| ) after the identifications <-> 2irp, dj. <-> 2wr^, and 



g f (k,k') =(p.(k-k') 



(35) 



Ref. 14 used the same fermionic state counting in the 
TBA as above but the bosonic statistical interaction 
function. With this remark, the identification ( |35| ) agrees 
with that obtained in Ref. 14 if one takes into account 
the relation (|27]) . But we stress that interpretation of the 
formula of reduction of the space of single-particle states 
( pi| ) in the TBA scheme of Refs. 25 and 26 necessar- 
ily implies that single-particle states are to be specified 
as fermionic. On the other hand, it is natural to expect 
that in an alternative, bosonic state counting in the TBA, 
which is discussed in Ref. 28, the reduction formula ( |2^ ) 
with the bosonic single-particle states should be used. 

As the first example, we consider the system of 
bosons with the two-body (5-function potential V(x) = 
c8{x). That this system models a possible statistics 
for identical jaarticles in one dimension was pointed, 
out long agoEj. The two-particle phase shift readsE3 
9(k) = -2tan- 1 (fc/c). Eq. (§5|) then yields g f (k,k') = 
—2c/L[c 2 + (k — k') 2 ]. In the limits c — > and c — > oo 
we have g f (k,k') = ^S(k - k') and g f (k,k') = 0, re- 
covering the thermodynamics for free bosons and free 
fermions as it should beEa. Note that after the identi- 
fications rife <-> 2irp and (|35|), Eq. ( |3^ ) coincides with 
the T = equation for the momentum, distribution for 
bosons with the (5-function interactior£3. 

Consider now the Calogero-Sutherland system ( |l9| ) 
(without the harmonic potential) . The Schrodinger equa- 
tion for the relative problem of two particles (—d 2 + \(\ — 
l)/.x 2 )-0 = k 2 tp, where x — x 2 — X\ and k = \{k x — k 2 ) 
are the relative coordinate and momentum, reduces to a 
Bessel equation so that any solution of the Schrodinger 
equation can be represented as a linear combination 
of the two solutions, ^/xJ\^i/2(kx) and yfxJi/2-\(kx), 
if A — i is a noninteger number (for A — h integer, 
the two above solutions are linear dependent, and an- 
other fundamental system of solutions, \/xJ\-i/2(kx) 
and \/xY\^i/2{kx), involving Bessel functions of the first 
and second kinds, should be used). These two solutions, 
for A > and A < 1, respectively, correspond to the 
solutions of classes I and II for the Calogero-Sutherland 
model discussed in the previous section. 



Consider the solution ^fxJ\-\/2{kx) for A > 0. For 
kx 3> 1, it has the asymptotics 7r -1 / 2 cos(fcx — J^A). 
Comparing this with the asymptotic wave function (|33|) 



for two particles, $ oc e 



S{k)e lkx , we get S(k) 



e~ for k > 0. In order to determine the phase shift 
9(k) = i\n(—S(k)) for all the values of A simultaneously, 
we consider the In z on its Riemann surface glued of the 
complex planes (k = 0,±1,±2,...) which are cut 
along their positive semiaxes and are characterized by 
2kn < argz < 2{k + 1)tt; in the D , the branch of the 
logarithm lnz = ln|z| + iargz with < argz < 2n is 
chosen. Then we obtain 8(k) = ir(\ — 1) for k > 0, and 
accounting for the antisymmetry of the phase shift, fi- 
nally, 6(k) = 7r(A — l)sgn(fc). This recovers the result 
obtained by SutherlandEj, our choice of the unique solu- 
tion of the scattering problem thus corresponds to that 
of Ref. 26. 

Eq. @ then yields g f (k,k') = 2f (X B - l)S(k - k') 
According to (p7|), the same statistics may be described 
by the bosonic statistical interaction function g b (k, k') = 
^X B 5{k~k'). 

Similar considerations show that the fcrmion Calogero- 
Sutherland system is equivalent to a free system with 



gt{k,k') 
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XpSik — k'). These results agree with 



those obtained in the previous section for the Calogero- 
Sutherland system in a harmonic well. 

In conclusion, we have proposed a formulation of the 
Haldane approach to statistics for spinless identical parti- 
cles considering a variable number of bosonic (fermionic) 
single-particle states which depends on the number of 
particles in the system. The variation of the num- 
ber of single-particle states is governed by the bosonic 
(fermionic) statistical interaction function which reduces 
to a single parameter for exclusion statistics [bosonic 
(fermionic) ES parameter]. Thermodynamic quantities 
for the ideal gas are determined by the momentum dis- 
tributions for single-particle states together with distri- 
butions of particles over single-particle states. The equa- 
tions for these two distributions in the most generic case 
are Eqs. (|29|)-(|30|). The fermionic state picture seems to 
be the most convenient for studying the low-temperature 
properties. 

The picture of variable number of single-particle state 
exactly corresponds to the description of states in the 
TBA for spinless particles. We also note that the pro- 
posed formulation of ES allows one to avoid encounter- 
ing negative probabilities which arise in the formulation 
of ES within the approach with nonvariable number of 
single-particle statestilc 2 ]. 

Generalization to several species of particles is possible 
and will be reported elsewhere. We also expect that the 
above formulation of exclusion and more general statis- 
tics can be extended to spinning particles where many- 
particle states of more general permutation symmetry 
than totally symmetric or antisymmetric ones should be 
used. 
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